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Abstract 

In this paper, we consider a planar dynamical system with a piecewise linear function 
containing an arbitrary number (but finite) of dropping sections and approximating 
some continuous nonlinear function. Studying all possible local and global bifur- 
cations of its limit cycles, we prove that such a piecewise linear dynamical system 
with k dropping sections and 2k + 1 singular points can have at most k + 2 limit 
cycles, k + 1 of which surround the foci one by one and the last, (k + 2)-th, limit 
cycle surrounds all of the singular points of this system. 
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1 Introduction 



The paper is based on the applications of Bifurcation Theory developed by An- 
dronov, Arnold, Thorn, Whitney, Zeeman et al. and can be used for modeling 
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problems, where system parameters play a certain role in various bifurcations. 
The theoretical studies of bifurcations deal with so-called universal problems. 
This means that sufficiently many parameters are available for universality of 
generic families of dynamical systems in the context at hand, under a relevant 
equivalence relation. This has led to the classification of generic, local bifur- 
cations. In many applications, models have a given number of parameters. 
Moreover, the bifurcation analysis, taking place in the product of phase space 
and parameter space, is not restricted to local features only. On the contrary, 
often the interest is the global organization of the parameter space regarding 
bifurcations which can be both local and global. 

This paper deals with so-called sewed dynamical systems, i.e., with systems 
for which the domain of definition is divided into sub-domains where different 
analytical systems are defined. The trajectories of these partial systems are 
sewed in one way or another on the boundaries of the sub-domains. Such sys- 
tems have some typical features, namely: 1) the system sewing is immediate 
from the physical meaning of the problem under consideration; 2) the system is 
piecewise linear, i.e., the partial systems from which it is sewed are linear sys- 
tems; 3) on the line of sewing, a point map (a first return function) is defined, 
what allows to determine the character of the system under consideration. 

Piecewise linear dynamical systems always contain some parameters and, un- 
der the variation of these parameters, the qualitative behavior of the systems 
can obviously change. We will consider the simplest bifurcations possible in 
the sewed systems when the sewing lines are unchanged under the parameter 
variations. It is natural to consider the following bifurcations which are similar 
to the simplest bifurcations of continuous dynamical systems: 1) the bifurca- 
tion of a singular point of focus type; 2) the bifurcation of an immovable point 
of focus type, a quasi-focus; 3) the bifurcation of a sewed limit cycle; 4) the 
bifurcation of a sewed separatrix going from a saddle to another saddle (the 
saddles can be both sewed and unsewed); 5) the bifurcation of a separatrix of 
a saddle-shaped singular point (sewed or unsewed) going from a saddle-shaped 
singular point to another such point or to a saddle (sewed or unsewed); 6) the 
bifurcation of a sewed saddle-node; 7) the bifurcation of a sewed separatrix of a 
saddle-node (sewed or unsewed) going out of the saddle-node and going back 
to it. Besides, some specific bifurcations can occur in sewed systems. Since 
in such systems, for example, arches of attraction or repulsion composed of 
immovable points can be similar to singular points, some bifurcations which 
are similar to the generation of a limit cycle from a focus can occur in the 
corresponding constructions. 

Piecewise linear systems have many applications in science and engineering. 
Special cases of such systems provide mathematical models for mechanical sys- 
tems with Coulomb friction, for valve oscillators with a discontinuous char- 
acteristic, for direct control systems with a two-point relay mechanism, for 
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planar dynamical systems modeling neural activity, etc. Despite their sim- 
ple structure and relevance to the applications, there is, to the best of our 
knowledge, no complete study of their dynamical properties. In most exist- 
ing papers, which deal with planar dynamical systems with piecewise linear 
right-hand sides, either the systems considered are continuous or only par- 
ticular cases are investigated. The first analytical results on such systems go 
back to Andronov, Vitt, and Khaikin in the 1930s (see p]). The existence 
and non-existence of an asymptotically stable periodic solution (limit cycle) 
of a piecewise linear system can be comparatively easily proved. However, for 
example, periodic solutions with sliding motion are of great importance to the 
applications. In particular, they describe the so-called stick-slip oscillations 
which appear in mechanical systems with dry friction. 

The main objective of the present paper is to provide a complete analysis of 
the dynamical properties of piecewise linear systems, their dependence on the 
system parameters studying, first of all, their limit cycle bifurcations. There 
are several ways to investigate the qualitative dynamics of such systems [7] . 
There are also numerous methods and good results on studying limit cycles. 
However, the most important impulse to their studying was given by the 
introduction of ideas coming from Bifurcation Theory [2], [3J, [SJ, [13] • We 
know three principal limit cycle bifurcations which we will study [8J: 1) the 
Andronov- Hopf bifurcation (from a singular point of centre or focus type); 

2) the separatrix cycle bifurcation (from a homoclinic or heteroclinic orbit); 

3) the multiple limit cycle bifurcation. All of these local bifurcations will be 
globally connected and will be applied to the qualitative analysis of piecewise 
linear dynamical systems. 



2 Preliminaries 

In this paper, geometric aspects of Bifurcation Theory are used and developed. 
It gives a global approach to the qualitative analysis and helps to combine 
all other approaches, their methods and results. First of all, the two-isocline 
method which was developed by Erugin is used [8] . An isocline portrait is the 
most natural construction for a polynomial equation. It is sufficient to have 
only two isoclines (of zero and infinity) to obtain principal information on the 
original polynomial system, because these two isoclines are right-hand sides of 
the system. Geometric properties of isoclines (conies, cubics, quadrics, etc.) are 
well-known, and all isoclines portraits can be easily constructed. By means of 
them, all topologically different qualitative pictures of integral curves to within 
a number of limit cycles and distinguishing center and focus can be obtained. 
Thus, it is possible to carry out a rough topological classification of the phase 
portraits for the polynomial dynamical systems and for the corresponding 
piecewise linear systems. It is the first application of Erugin's method. After 
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studying contact and rotation properties of isoclines, the simplest (canonical) 
systems containing limit cycles can be also constructed. Two groups of the 
parameters can be distinguished in such systems: static and dynamic. Static 
parameters determine the behavior of phase trajectories in principle, since they 
control the number, position and character of singular points in a finite part of 
the plane (finite singularities). The parameters from the first group determine 
also a possible behavior of separatrices and singular points at infinity (infinite 
singularities) under the variation of the parameters from the second group. The 
dynamic parameters are rotation parameters. They do not change the number, 
position and index of the finite singularities, but only involve the vector field in 
the directional rotation. The rotation parameters allow to control the infinite 
singularities, the behavior of limit cycles and separatrices. The cyclicity of 
singular points and separatrix cycles, the behavior of semi-stable and other 
multiple limit cycles are controlled by these parameters as well. Therefore, 
by means of the rotation parameters, it is possible to control all limit cycle 
bifurcations and to solve the most complicated problems of the qualitative 
theory of dynamical systems (both continuous and piecewise linear). 

In [S] some complete results on continuous quadratic systems have been pre- 
sented, and some preliminary results on generalizing geometric ideas and bi- 
furcation methods to cubic dynamical systems have been obtained. So, in [TT] . 
a canonical cubic dynamical system of Kukles-type was constructed and the 
global qualitative analysis of its special case corresponding to a generalized 
Lienard equation was given. In particular, it was proved that the foci of such 
a Lienard system could be at most of second order and that the system could 
have at least three limit cycles in the whole phase plane. Moreover, unlike all 
previous works on the Kukles-type systems, by means of arbitrary (including 
as large as possible) field rotation parameters of the canonical system, the 
global bifurcations of its limit and separatrix cycles were studied. As a re- 
sult, the classification of all possible types of separatrix cycles was obtained 
and all possible distributions of limit cycles were found for the generalized 
Lienard system. In [5], the global qualitative analysis of centrally symmetric 
cubic systems which are used as learning models of planar neural networks 
was established. All of these results can be generalized to the corresponding 
piecewise linear Lienard-type dynamical systems. 

Such systems have been considered in several paper. For example, Bautin [3] 
studied a dynamical system which was used in radio-engineering for describing 
tunnel diode circuits, where the nonlinear function was approximated by a 
piecewise linear function composed of three linear pieces. Giannakopoulos and 
Pliete [12] generalized that results for a piecewise linear system with a line 
of discontinuity modeling stick-slip oscillations which appeared in mechanical 
systems with dry friction. In this paper, the obtained results are generalized 
to an arbitrary piecewise linear Lienard-type system, where the polynomial 
is approximated by a piecewise linear function composed of an arbitrary (but 
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finite) number of linear pieces. In particular, the following bifurcations are 
studied for such a system: 1) the bifurcation of singular points (division of the 
parameter space according to the number and character of singular points, 
stability of singular points lying on the lines of sewing, generation of limit 
cycles from singular points of focus type under transferring the points through 
the lines of sewing, generation of limit cycles (hyperbolic and semi-stable) from 
the boundaries of the domains filled by closed trajectories); 2) the separatrix 
bifurcation (location of the bifurcation curves for the separatrix loops, stability 
of the separatrix loops); 3) the bifurcation of multiple limit cycles (location of 
the bifurcation curves for the limit cycles of various multiplicity, qualitative 
structure of the division in the phase plane). 



3 Main Results 

Consider the system 

x — y — (p(x), y = {3 — ax — y, a > 0, (3 > 0, (1) 

where (f(x) is a piecewise linear function containing k dropping sections and 
approximating some continuous nonlinear function. The line (3 — ax — y = 
and the curve y = <p(x) can be considered as the isoclines of zero and infinity, 
respectively, for the corresponding equation. Such systems and equations may 
occur, for example, when tunnel diode circuits and some other problems are 
studied (see [I], [3], @], [7], [13]). 

Suppose that the ascending sections of (pQ) have an inclination ki > and the 
descending (dropping) sections have an inclination k 2 < 0. Then the phase 
plane of ([1]) can be divided onto 2k + 1 parts in every of which is a linear 
system: the ascending sections are in k + 1 strip regions (/, V, . . . , 2K + 1) 
and the descending sections are in other k such regions (II, IV, VI, . . . , 2K). 
The parameters k%, k 2 , and also a can be considered as rotation parameters 
for the sewed vector field of ([T]) (see [I], [S]). 

System (pQ) can have an odd number of rough singular points: 1, 3, 5, ... , 2k + l. 
If ((I]) has the only singular point, this point will be always an antisaddle 
(center, focus or node). A focus (node) will be always stable in odd regions 
and unstable in even regions ii k 2 > 1. If system ([1]) has 2k + 1 singulari- 
ties, then k of them are saddles (they are in even regions) and k + 1 others 
are antisaddles (foci or nodes) which are always stable (they are in odd re- 
gions). The pieces of the straight lines (3 = x 2 %-\a + yn-\ and (3 = x^a + y 2 % 
[i = 1,2, ... ,k), where (a^t-i, 2/2i-i) and [x 2 i,y2i) are the coordinates of the 
upper and lower corner points of the curve (f(x), respectively, form a discrim- 
inant curve separating the domains in the plane (a, (3), where a < k 2 , with 
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different numbers of singular points. The points of the discriminant curve cor- 
respond to the sewed singularities of saddle-focus or saddle-node type (a < k 2 ) 
and its corner points correspond to the unstable equilibrium segments (a = k 2 ) 
which coincide with the dropping sections of the curve y = (p(x). In the case 
when k 2 < 1, closed trajectories cannot exist and only bifurcations of singular 
points are possible in system (TjQ). Therefore, we will consider further only the 
case when k 2 > 1 and {k\ — l) 2 < 4k 2 giving various bifurcations and, first of 
all, the bifurcations of limit cycles. 

In Section 4, studying local bifurcations of limit cycles, we prove that: 1) the 
only limit cycle can appear in system ([1]) under a displacement of a singular 
point of focus type through the sewing line of (Q, and this limit cycle will be 
unique in the corresponding pair of the strip regions; 2) at most two limit cycles 
can be generated by the boundary of the domain filled by closed trajectories 
of system ([I]), and these limit cycles can be only outside the boundary; 3) the 
equilibrium segment (the dropping section) of ([1]) can generate at most two 
unstable limit cycles surrounging two stable foci one by one which appear 
simultaneously from the ends of the segment; 4) a nondegenerate separatrix 
cycle of system (pQ) with 2k + 1 simple singular points can generate at most 
k + 1 unstable limit cycles inside its loops (digons) or the only unstable limit 
cycle outside it. In Section 5, using the obtained results and applying the 
Wintner-Perko termination principle for sewed multiple limit cycles, we sudy 
global bifurcations of limit cycles and prove that system (TjQ) with k dropping 
sections and 2k + 1 singular points can have at most k + 2 limit cycles, k + 1 
of which surround the foci one by one and the last, (k + 2)-th, limit cycle 
surrounds all of the singular points of ([!]). 



4 Local Bifurcations of Limit Cycles 

Following [3J, we begin with studying local bifurcations of limit cycles. The 
limit cycle of system (TjQ) will be called small if it belongs to at most two 
adjoining regions; the cycle will be called big if it belongs to at least three 
adjoining regions. 

4-1 The Limit Cycle Bifurcation from a Singular Point of Focus Type under 
a Displacement of the Singular Point through the Line of Sewing 

Lemma 4.1. Under a displacement of a singular point of focus type through 
the line of sewing of system at most one limit cycle can appear in this 
system, and this limit cycle will be unique in the corresponding pair of the strip 
regions of £!]]■ 
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Proof. Let us study first stability of a singular point on the line of sewing. 
Suppose that the straight line j3 — ax — y = passes through the corner point 
(xi,yi) of the curve y = ip(x) on the boundary of regions I, II and that 
a > (k 2 + l) 2 /4. Then region / (II) will be filled by pieces of trajectories of 
a stable (unstable) focus of system ([I]). 

Introduce positive coordinates So (lower (xx,yi)) and Si (upper (xi,yi)) on 
the line of sewing of regions / and II; S 2 (lower (x 2 , 2/2)) and S3 (upper (x 2 , 1/2)) 
on the line of sewing of regions II and III, etc. The maps Sq — > S± along the 
trajectories of region / and S\ — > Sq along the trajectories of region 77 are 
written as follows: 

S l = S e 1U71,ul \ Sq = Ste™ 2 ^ 2 , (2) 

where <7j, u>i (i = 1,2) are the real and imaginary parts of the roots of the 
characteristic equation for a singular point of regions /, /J, respectively. 

The singular point (x\ , y\) will be a sewed center (Sq = So) iff o\ / u\+a2 /uj 2 = 0, 
i.e., when a = a* = (1 — k\/kz)/(k2 — k\ + 2). The sewed focus (x\,yi) will 
be stable (So < So) when a > a* and unstable (So > So) when a < a*. 

Consider the return map Sq — > Sq along the trajectories of regions / and //. 
For region J, we will have 

50 = (uji coscjiri — (Ji sincJiTi — Ldie~ aiTl ) = 5o((t\), 

51 = (uJi cosujiTi + <7i sintUxTx — uJie aiTl ) = 5ox( r i), 



where 5q is the distance from the boundary of regions /, // to the singular 
point; ( and x are monotonic functions. The return map along the trajectories 
of region II has a similar form. 



Calculation of the first derivative for the return map gives 

e 2(<71Tl+CT2T 2 ) /^\ 



dS _ ^0 ^2(a 1 r 1 + CT2 r 2 ) 



dS S 



where (i = 1,2) is motion time along the trajectories of regions J, /J, 
respectively; cij = (1 + k^/2 (i = 1,2). 

Suppose that the singular point of system ([1]) is in region /. Then, along a 
periodic solution of Q (when S = S ) and for increasing S , time T\ decreases 
(to the value n/ui), time r 2 increases (to the value tt/^), and derivative 
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fll]) also increases. Therefore, at most two intersection points of the graph of 
the return map S — > S with the bisectrix of the plane (S ,S ) can exist; 
moreover, the fixed point with the smaller coordinate must be stable and the 
fixed point with the grater coordinate must be unstable. Since, under our 
assummption, the singular point is in region / and it is a stable focus which 
cannot be surrounded by a stable limit cycle, then at most one limit cycle can 
exist in regions / and //; moreover, this limit cycle can be only unstable. 

Suppose that the singular point of system ([1]) is in region II. Then, for in- 
creasing 5*o, time T\ increases and time r 2 decreases. In a similar way, we can 
show that at most one (stable) limit cycle can exist in this well. 

Let the straight line f3 — ax — y = pass through the upper corner point 
(xi,yi) of the curve y = <p(x) again and suppose that a > a* (the case when 
a < a* is considered absolutely similarly). The sewed focus is stable in this 
case. The trajectory passing through the lower corner point (£2,2/2), by ([2]), 
tends to the singular point for t — > +00. This trajectory remains a spiral under 
a small displacement of the line (3 — ax — y = 0. If, under such a displacement, 
the singular point gets into region II, then it becomes unstable and, hence, 
at least one stable limit cycle appears in system ([1]). As follows from above, 
this cycle is unique. If, under a displacement, the singular point gets into 
region J, obviously, limit cycles cannot appear in system ([1]). The same results 
can be obtained for regions 77/ and IV, . . . , 2K — 1 and 2K. The lemma is 
proved. □ 



4-2 The Limit Cycle Bifurcation from the Boundary of the Domain Filled by 
Closed Trajectories 



Lemma 4.2. The boundary of the domain filled by closed trajectories of system 
(T7J) can generate at most two limit cycles, and these limit cycles can be only 
outside the boundary. 



Proof. Consider now the map S = /(So) sewed of two pieces: S = £(S ) 
along the trajectories in regions I, II, . . . , 2K and S = ^(So) along the 
trajectories in all regions, I, II, . . . , 2K, 2K + 1. The map So — > Si in region 
I is given by (j3J). The maps Si -> S 3 , S 3 — ► S 5 , . . . , S 2 k-i -> S 2 k-2 (S 2 k-i -> 
S2fc+i> S2fc+i — > S2fc, S2A; — > S2k-2), S 2 k-2 —* S 2 k-4, ■ ■ ■ , S2 — > So have similar 
forms. 



The derivatives for the functions £(So), tp(So) are given by the following ex- 
pressions, respectively: 



_ £0 ^ ai (j 1+ T++T-+...+T 2k _ 1 )+a 2 {T++T-+... + T+ k _ 2 + 

dS S 



2k-2^'2k-2' 



■ 



(5) 
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dS, 



So 
So 



e 2( CT1 ( T1 +r3++r3 +... +T2fe+1 ) +0 - 2 ( r ++T 2 +...+T+ +r 2fe )) 



(6) 



where 7i, r 2 fc_i, r 2 fc + i are motion times in regions J, 2if — 1, + 1 and 
( r 2i); r 2i+i ( r 2i+i)> * = 1) 2, • • • , k, are motion times in the upper (lower) parts 
of regions II, III, . . . , 2K, respectively. 

Let Sq = Sq be the boundary value separating the intervals of definition of the 
maps £(So) an d if>(So). By ((SJ) and (Q, the function S = /(So) is differentiable 
at the point Sq. 

Suppose now that the straight line (3 — ax — y = passes through the corner 
point (xi,yi) and that a = a*. Let us show that system has no limit cycle 
in this case. 

The return map in the plane (So, So) is sewed of two pieces: the bisectrix seg- 
ment So = So < Sq and the curve So = if>(So). Since the function 
So — f(So) is differentiable at the sewing point Sq, then dSo/dS = 1 for 
a = a * (besides, if follows from ® that d 2 ^ /d^ < 0). For increasing S 
(from the value Sq), the exponent degree in ([6]) will monotonically decrease at 
the sewing point (from the zero value). It is clear that the curve So = if)(So) 
cannot have other points of intersection (or tangency) with the bisectrix in 
addition to So = Sq. For So > Sq, the curve So = ip(So) lies below the 
bisectrix. Therefore, the spirals sewed of the trajectories of ([1]) in regions 
/, 77, . . . , 2K, 2K + 1 tend to the boundary So in regions /, //, . . . , 2K. 

For a small variation of the parameters a and f3, the return map graph of the 
changed system will be in a small neighborhood of the return map graph of 
the original system. If we move along the half-line L\ = (L\ = (3 — ax\ — y^ 
a > k,2) from the value a = a* to the side of decreasing a, then the return map 
will be either the straight line passing through the origin above the bisectrix 
(for S < Sq) or the curve S = if) (So) intersecting the bisectrix one time (for 
So > Sq). In this case, the only stable limit cycle occurs from the boundary of 
the domain filled by closed trajectories. For further decreasing (3, the return 
map graph will move from the origin along the So-axis and will intersect the 
bisectrix two times (an unstable limit cycle occurs from the focus under its 
displacement through the line of sewing). If we move along the half-line to 
the side of increasing a from the value a = a*, decreasing (3 after that, then 
the return map graph will be completely below the bisectrix. It follows from 
continuity and differentiability of the return map that there exists a point 
(a*, (3*) on the half-line L\ = at which a bifurcation curve of double limit 
cycles is created and then extended into some half-neighborhood of this point 
(below the half-line) bounding a domain with two limit cycles in the parameter 
space. 
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The points of this bifurcation curve correspond to the tangency points of 
the return map graph with the bisectrix of the plane (So, S ). Such tangency 
is impossible for S < Sq. Therefore, a double limit cycle can appear only 
from the boundary of the domain filled by closed trajectories of (pQ) when 
So — Sq. Thus, we have proved that at most two limit cycles can be generated 
by the boundary and that these two limit cycles can be only outside the 
boundary. □ 

4-3 The Limit Cycle Bifurcation from the Ends of the Equilibrium Segment 

Lemma 4.3. The equilibrium segment (the dropping section) of system (QJ) 
can generate at most two unstable limit cycles surrounging two stable foci one 
by one which appear simultaneously from the ends of the segment. 

Proof. Suppose that a part of the straight line j3 — ax — y = coincides 
with a dropping section of (jTJ), for example, with the first one (a = k 2 ). The 
dropping section of ([I]) will be an unstable equilibrium segment and regions 
/, II (because of the condition {k\ — l) 2 < 4A; 2 ) will be filled by trajectories 
of the stable foci. It is easy to obtain an explicit expression for the map of the 
half-line So into itself: 

S = S e 2 ™ 1 ^ 1 + 5{k 2 - 1)(1 + e lfflM ), (7) 

where 5 is the width of regions II. This map has the only stable fixed point. 

Rotate the straight line /3 — ax — y = counterclockwise around some point 
on the dropping section. The equilibrium segment will be destroyed and three 
singular points will appear: a saddle in region 77 and two stable foci in regions 
I and III. Let a = k 2 — s, where e is positive and sufficiently small. Being 
restricted by the first degree of e, we can obtain the following separatrix slopes: 
— l + e/(k 2 — 1) for the a-separatrix and —k 2 — e/(k 2 — l) for the a;-separatrix. 

For a = k 2 , the trajectories going out the point, where a saddle appears 
when e ^ 0, tend to a big stable limit cycle surrounding all of the singular 
points; for e > 0, the a-separatrices of the saddle from region II lie in a small 
neighborhood of the trajectories going out the same point for e = and, hence, 
they also tend to the same limit cycle. Therefore, the tu-separatrices can just 
go out unstable limit cycles lying in regions I — II, II — III and surrounding 
stable foci appearing in regions / and II under the rotation of the straight 
line j3 — ax — y = 0. 

Thus, under the rotation of this straight line, unstable limit cycles surrounding 
stable foci appear from the ends of the equilibrium segment (the cycles and the 
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foci appear simultaneously). There is the only limit cycle in a neighborhood 
of each of the foci. This follows from that the derivative constructed by using 
the trajectories of the saddle from region II will be also given by expression 
01]), with the only difference: for increasing So, r 2 — > +00. The same result can 
be obtained for all of the other dropping sections. This completes the proof of 
the lemma. □ 



4-4 The Separatrix Cycle Bifurcations 

The simplest type of separatrix cycles of p]) is a so-called eight-loop formed 
by two ordinary saddle loops. In the case of 2k + 1 simple singular points, a 
separatrix cycle can contain k + 1 saddle loops, the first and the last of which 
are ordinary loops with one rough saddle on each and the k — 1 others are 
separatrix digons with two rough saddles on each. Such a separatrix cycle will 
be called nondegenerate. In the cases when the straight line (3 — ax — y = 
passes through the corner points of the curve y = (p(x), we will have degenerate 
separatrix cycles of lips-type containing one or two sewed saddle-nodes. Let us 
consider the bifurcations of nondegenerate separatrix cycles. To consider the 
cases of degenerate separatrix cycles, we have either to generalize the results 
of [6] or to apply the Wintner-Perko termination principle [8] , [H] which will 
be done in the next section. 

Lemma 4.4. A nondegenerate separatrix cycle of system (TJP with 2k+l simple 
singular points can generate at most k + 1 small unstable limit cycles inside 
its loops (digons) or the only big unstable limit cycle outside it. 

Proof. Let us show first that the bifurcations of a separatrix loop do not 
depend on the parameter (3. Fix the parameters fci, k 2 , a and suppose that 
the line (3 — ax — y = passes through the upper corner point (xi,yi) of the 
curve y = <p(x) for some value j3 = (3q. Varying (3 by the value k (k = f3o — (3) , 
we will show that a separatrix loop cannot be formed under such a variation. 
Let S' and S[ be the segments intercepted on the boundary of regions I and 
77; So and Si be the coordinates for return map ([3]) on the same boundary. It 
follows from that 



where ( 1 is the inverse function for (. The values o~\ and uj\ do not depend 
on (3\ therefore, the functions x an d C a l so do not depend on it. 

Since y = <f(x) is a piesewise linear function, then, under the variation of (3, 



Si = S oX C\S /S ) 




11 



the values S' , S[, and 5o will be proportional to k 



(9) 



Sewing the trajectories on the boundary of regions I and II (assuming S'q = S ), 
we can find from (jSJ), (jUJ): 

Si = 7i^xC~ 1 (7o/7i) = 73^, (10) 

Si/S[ = 73/72 = const. (11) 



Thus, under the fixed ki, k 2 , and a, the values Si and S[ have a constant 
ratio; hence, a separatrix loop (Si = S[) cannot be formed under a variation 
of the parameter (3. The same conclusion can be made for any other type of 
the separatrix cycles of (jTj). It means that separatrix cycles can be formed or 
destroyed only under a variation of the parameter a for fixed k\ and hi. 

As is known [I], the character of stability of a sewed nondegenerate separa- 
trix cycle is determined by the sign of its saddle quantities which are always 
positive in our case, when the saddles are inside the boundary of even regions 
II, IV, . . . , 2K and k 2 > 1 (Theorems 44 and 47 from [2] are valid for the 
piecewise linear dynamical systems as well). It follows that a nondegenerate 
separatrix cycle of (pQ) is always unstable (inside and outside) and, under a 
variation of a, such a separatrix cycle can generate at most k + 1 small un- 
stable limit cycles inside its loops (digons) or the only big unstable limit cycle 
outside it when system ([1]) has 2k + 1 simple singular points. The lemma is 
proved. □ 



5 Global Analysis 



For the global analysis of limit cycle bifurcations, we will use the Wintner- 
Perko termination principle which is formulated for planar relatively prime 
analytic systems connecting the main local bifurcations of limit cycles j8], |14j . 



Theorem 5.1 (Wintner— Perko termination principle). Any one-para- 
meter family of multiplicity-m limit cycles of a relatively prime analytic system 

x = f{x, 11), (12) 



where x 6 R 2 , /x e R n , f e R 2 (/ is an analytic vector function ), can be 
extended in a unique way to a maximal one-parameter family of multiplicity-m 
limit cycles of ( TUf) which is either open or cyclic. 
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// it is open, then it terminates either as the parameter or the limit cycles 
become unbounded; or, the family terminates either at a singular point of HUty , 
which is typically a fine focus of multiplicity m, or on a (compound) separatrix 
cycle of (flDj). which is also typically of multiplicity m. 

As was shown in jl], the return map constructed in a neighborhood of a 
multiple limit cycle of a piecewise linear dynamical system (under some general 
assumptions which are valid in our case) will be an analytic function; therefore, 
all statements formulated for multiplicity-m limit cycles of analytic system 
(1121) will be also valid for sewed multiplicity-m limit cycles. Applying the 
Wintner-Perko termination principle and the results on local bifurcations of 
limit cycles obtained in the previous section for system (pQ), we will prove the 
following theorem. 

Theorem 5.2. System (QP with k dropping sections and2k + l singular points 
can have at most k + 2 limit cycles, k + 1 of which surround the foci one by 
one and the last, (k + 2)-th, limit cycle surrounds all of the singular points 

of [ID- 

Proof. Let us consider global bifurcations of limit cycles of p]). Suppose again 
that the zero isocline (3 — ax — y = passes through the corner point (x±, yi) of 
the infinite isocline y = (f(x) and that a > a*. In this case, the only singular 
point in the phase plane is a sewed stable focus and all trajectories of (pQ) tend 
to it when t — > +oo. For decreasing a {k 2 < a < a*), the sewed focus becomes 
unstable and a stable limit cycle is generated from the boundary curve of the 
domain with closed trajectories (immediately after passing the value a* by 
the parameter a). 

For a = k 2 , the first dropping section of ([I]) will coincide with a part of the 
straight line (3 — ax — y = and an unstable equilibrium segment will appear 
inside the stable limit cycle. If we rotate the line (3 — ax — y = around an 
interior point of the segment (changing both of the parameters, a and 
two unstable limit cycles surrounding stable foci one by one will be generated 
from the ends (x\,yi) and (x 2 ,y 2 ) of the equilibrium segment. Under the fur- 
ther rotation of the line f3 — ax — y = 0, it will pass first through the next 
corner point, (24, 2/4), and then, successively, through the points (xe,ye), . . . , 
(x2k,y2k)- Every time, the corner point becomes a sewed saddle-node gener- 
ating an unstable limit cycle surrounding a stable focus. So, we will get a 
piecewise linear system with 2k + 1 singular points having at least k + 1 small 
unstable limit cycles surrounding the stable foci one by one inside a big stable 
limit cycle, k + 2, surrounding all of the singular points. 

Under the further rotation of the zero isocline, all k + 1 small limit cycles 
simultaneously disappear in a nondegenerate separatrix cycle consisting of 
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k + 1 loops (digons), this separatrix cycle generates a big (unstable) limit cycle 
which combines with another big (stable) limit cycle of ([1]) forming a semi- 
stable (double) limit cycle which finally disappears in a so-called trajectory 
condensation. 

Let us prove now that system ([1]) cannot have more than k + 2 limit cycles. 
The proof is carried out by contradiction. Since a small limit cycle is always 
unique in the corresponding strip regions (Lemma 4.1), suppose that system 
(CQ) with three field rotation parameters, k\, k 2 , and a, has three big limit 
cycles. Then we get into some domain in the space of these parameters which 
is bounded by two fold bifurcation surfaces forming a cusp bifurcation surface 
of multiplicity-three limit cycles [8], [14"] . 

The corresponding maximal one-parameter family of multiplicity-three limit 
cycles cannot be cyclic, otherwise there will be at least one point corresponding 
to the limit cycle of multiplicity four (or even higher) in the parameter space. 
Extending the bifurcation curve of multiplicity-four limit cycles through this 
point and parameterizing the corresponding maximal one-parameter family of 
multiplicity-four limit cycles by a field rotation parameter, for example, by 
the parameter a, we will obtain a monotonic curve (see [8], [2]) which, by 
the Wintner-Perko termination principle (Theorem 5.1), terminates either at 
the boundary curve of the domain filled by closed trajectories of (CQ) or on 
some degenerate separatrix cycle of ([I]). Since we know at least the cyclicity 
of the boundary curve which is equal to two (Lemma 4.2), we have got a 
contradiction with the termination principle (Theorem 5.1) stating that the 
multiplicity of limit cycles cannot be higher than the multiplicity (cyclicity) 
of the end bifurcation points in which they terminate. 

If the maximal one-parameter family of multiplicity-three limit cycles is not 
cyclic, using the same principle (Theorem 5.1), this again contradicts with 
the cyclicity result for the boundary curve (Lemma 4.2) not admitting the 
multiplicity of limit cycles higher than two. Moreover, it also follows from the 
termination principle that the degenerate separatrix cycles of (CQ) cannot have 
the multiplicity (cyclicity) to be higher than two. Therefore, according to the 
same principle, there are no more than two big limit cycles in the exterior 
domain outside the boundary curve of ([Tj). 

The same results can be obtained by means of the new geometric methods 
developed in [9], [TO]. The phase portraits and bifurcation diagrams for system 
(CQ) will be similar to that which were constructed in [3J, [I]. Thus, system (TJQ) 
with 2k + 1 singular points cannot have more than k + 2 limit cycles, i.e., 
k + 2 is the maximum number of limit cycles of such system and the obtained 
distribution (k+1 small limit cycles plus a big limit cycle) is the only possibility 
for their distribution. This completes the proof of the theorem. □ 
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